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Abstract. For a class of random matrix ensembles with corre- 
lated matrix elements, it is shown that the density of states is 
given by the Wigner semi-circle law. This is applied to effective 
Hamiltonians related to the Anderson model in dimensions greater 
than or equal to two. 



It is a classical theorem due to Wigner [Wig] that the density of 
states of a growing sequence of real symmetric matrices with inde- 
pendent entries converges in distribution to the semi-circle law. More 
precisely, this means the following: Consider an n x n random matrix 
= (^7= o, n (p, q))i<p,q< n where, apart from the symmetry condition, 
the entries a n (p,q) are independent centered random variables with 
unit variance (and a growth condition on their moments). Then the 
expectation value of the moment of X n of order k > satisfies 



The coefficients are called the Catalan numbers, and the terms 
on the r.h.s. of (1) are precisely the moments of the semi-circle law, 
namely the absolutely continuous measure supported by the interval 
[-2,2] with density - x 1 . 

In this note we address the question: How many independent ma- 
trix elements are needed such that the above asymptotic behavior still 
holds? In other words, how many correlations may the matrix ele- 
ments have? This is similar to an extension of the classical central 
limit theorem to sums of random variables with decaying correlations. 
We consider here a particular class of such correlations, including cases 
with supplementary symmetries of the matrix elements (apart from 
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a n{p,<l) = a n(q,p))- This is motivated by applications, in particular a 
simplified effective Anderson model. For the two-dimensional effective 
Hamiltonian, the so-called "flip- matrix" model introduced in [Poi], it 
was shown in [BMR] that the density of states is semi-circular using 
super symmetric functional integrals. Bellissard has conjectured [Bel] 
that this result also holds for the effective Hamiltonians in higher di- 
mensions. As a consequence of our result we obtain a proof of this 
conjecture (see Section 2 for details). In Section 3, we indicate that 
known results on asymptotic freeness [Voi, Spe] extend to the present 
ensembles. This allows to calculate higher correlation functions. 

The precise set-up we consider is as follows. For each n G N, sup- 
pose given an equivalence relation ~ n on pairs P = (p, q) of indices in 
{1 , . . . , n} x2 . The entries of the matrix X n = (-^= a n (p, q))i< p , q < n are 
complex random variables, with a n (pi, qi), . . . , a n (pj, qj) independent 
whenever (pi, qi), . . . , (pj, qj) belong to j distinct equivalence classes of 
the relation ~ n . Just as above, the entries are supposed to be centered, 
to have unit variance, and to obey a moment condition, namely 

(2) m k : = sup max E(\a n (p, q)\ k ) < oo , 

n P,<2=1,---,™ 

for all k G N. For equivalent pairs (p,q) ~„ {p',q'), the relation be- 
tween a n (p, q) and a n (p', q') is not specified, and these variables may be 
correlated. For instance, we might have a n (p, q) = a n (p', q'), describing 
an auxiliary symmetry as in [Poi, BMR, Bel], but this is not necessary. 
However, the following conditions are imposed on the equivalence rela- 
tion: 

(CI) max#{( (Z ,p',( Z ')G{l,...,n} xS |( P ,(z)~ n (p',9')} = o(n 2 ) 
v 

(C2) max#{g , e{l,...,n}|(p,g)~ n (p',g')} < B 

p,q,p 

for some constant B < oo, and 

(C3) #{(p,g,p / )e{l,...,n} x3 |(p,g)~ r^ (g,p / )&P^P , } = o(n 2 ) . 
Finally, we require that 

(3) a n (p,q) = a n (q,p) , 

so that X n is Hermitian. Thus, (p,q) ~ n (q,p) for every p, q, which is 
consistent with (CI) to (C3). 

Theorem 1. The density of states of an ensemble of Hermitian random 
matrices X n obeying the moment bound (2) and associated to equiva- 
lence relations ~ n satisfying (CI), (C2) and (C3) is a semicircle law, 
that is, (1) holds. 
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Remark: The first condition (CI) means that the size of the equiva- 
lence classes, i.e., the number of entries correlated with a given entry 
of the matrix, cannot grow too fast. The second condition (C2) is more 
technical, but plays a key role in the proof. It could be relaxed, re- 
placing B by 0{n £ ) for every e > 0, provided (CI) were strengthened, 
replacing o(n 2 ) by o(n 2_£ ) for all e > 0. The third condition (C3) 
ensures, among other things, that there are not too many dependent 
rows, which could lead to many eigenvectors with eigenvalue zero. (An 
example illustrating this phenomenon is described below.) 

Our proof of this theorem is an extension of Wigner's original proof 
for the fully independent case [Wig, Spe, HP], the difference being that 
we use an equivalence relation on the pairs of indices rather than on 
the indices themselves (which may make Wigner's proof even more 
transparent to some readers). 

Proof of Theorem 1. To begin, let us call a sequence (Pi, . . . , P k ) 
of pairs "consistent" if neighboring pairs Pi = (pi,qi) and P; + i = 
(pi+i, qi+i) satisfy q\ = pi+i, where I = 1, . . . , k and k + 1 is cyclically 
identified with 1. Then, for any given k G N, we have 

(4) E - Tr„(X*) = — T Yl E MPi) • • • a n (P k )) , 

^ ft,... A =(i,i) 

where the sum runs over consistent sequences only. The terms of this 
sum are uniformly bounded as n — > oo: 

|E(a n (Pi) . . .a n (P k ))\ < m k , 

as follows from (2) and the Holder inequality. 

To group similar terms, let us associate a partition ir of {1, . . . , k} 
to each sequence (Pi, . . . , P k ) by means of 

/ ~ w m Pi ~ n P m , 

and refer to (Pi, . . . , P k ) as a "tt consistent sequence". Introducing the 
notation S n (-n) for the set of n consistent sequences with indices in 
{1, . . . , n}, equation (4) becomes 

(5) E^Tr n (X n fc ) 

= -HT E E E(a n (P 1 ).--a n {P k )), 

U TrGP(fc) (P 1 ,...,P k )ES n (n) 

where V(k) denotes the set of all partitions of {1, . . . , k}. 
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Let #7r be the number of blocks of the partition n. If #7r > |, 
then there has to be one singleton block, namely a block with only one 
element. Thus for any 7r consistent sequence the corresponding entry 
in a n (Pi) ■ ■ ■ a n (Pk) appears only once. As this variable is centered and 
independent of all others appearing, E(a n (Pi) • • -a n (Pk)) = 0. There- 
fore the sum in (5) can be restricted to partitions with less than or 
equal to | blocks. 

Next we argue that partitions with r = #7r < | give vanishing 
contribution in the limit n — > oo. Since the corresponding term in the 
sum over V(k) in (5) is bounded by x #S n (7r), it suffices to derive an 
upper bound on the number of n consistent sequences. Let us begin by 
choosing the pair P 1 . There are n 2 choices to be made, corresponding 
to n choices for each of the indices. Now, the first index of P 2 is fixed 
due to consistency. If 2 ~„- 1, then the second index can take at most 
B values due to condition (C2), otherwise it is unconstrained, and can 
take at most n values. Similarly, once we get to p, I < k, we have 
either at most B or at most n possible choices for the second index, 
depending on whether I ~ n j for one of j = 1, ...,/ — 1. For the last 
pair Pj, there is no freedom, due to consistency. This shows 

#S n (vr) < n 2 n r - 1 B k - r . 

Since r < |, this contribution is negligible compared to the prefactor 
l/n 1+k / 2 in (5) in the limit n — > oo. In particular, this implies that the 
limit in (1) vanishes if k is odd. 

For even k, we may now focus on the contributions coming from par- 
titions with exactly k/2 blocks. Moreover, if the partition has a single- 
ton the contribution vanishes as described above. Hence the partition 
has to be a pair partition, i.e., each block has exactly two elements. 
(Note that the "pair" in "pair partition" has nothing to do with the 
fact that each point I is associated to a pair p of indices.) 

Let VV(k) denote the pair partitions of k. To prove Theorem 1, we 
must control somewhat more carefully the asymptotics in n of #S n (n) 
for a pair partition n G W(k). This will be accomplished in the 
following four lemmas. 

Lemma 1. Let n G W(k) be a pair partition containing a pair of 
neighbors, that is m ~ n m + 1 for some m, and let ir' G W{k — 2) 
be the partition obtained by eliminating the corresponding pair (and 
relabeling Z i — > Z — 2 for m + 2 < I < k) . Then 

(6) #S' n (7r) < n x #S' n (7r / ) + o(nt +1 ) . 
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Proof. Let us look at the situation close to m and m+1. The indices are 
(p m -i,9m-i), (p m ,Qm), (Pm+i, Qm+i) and (p m+2 ,q m+ 2)- By consistency, 
we have q m = p m +i- Now consider separately the two cases (i) p m = 
q m+l and (ii) p m ^ q m+1 . 

In case (i), after two applications of the consistency condition, we get 
q m -i = Pm+2- Hence after eliminating the pairs P m and P m +i, we have 
a it' consistent sequence. Therefore, in this case there are n choices 
for q m = p m +i and at most #S n (ir') choices for the remaining indices, 
giving the first term in (6). 

In case (ii), by (C3) there are only o(n 2 ) choices for the triple p m , 
qm = Pm+i, qm+i- Since there are (k — 2)/2 pairs remaining, there 
are no more than n^ k ~ 2 ^ 2 B^ 2 ^ 2 ~ l choices for the remaining indices. 
To see this, start at m + 2 and proceed through the remaining indices 
cyclically as in the argument to eliminate partitions with #7r < k/2 
above. Combining the two factors o(n 2 ) x 0(n^ k ~ 2 ^ 2 ) gives the second 
term in (6). □ 

A partition 7r is called crossing if there are positions mi < m 2 < 
m 3 < m 4 such that mi ~ n m 3 and m 2 ~7r m 4 , otherwise it is called 
non-crossing. 

Lemma 2. Let n G W(k) be a crossing pair partition, then 

Proof. The first step is to apply Lemma 1 as many times as possible, 
eliminating all nearest neighbor pairs in n and the resulting reduced 
partitions. The end result is that 

#Sn(vr) < n r #S n (n') + o(n^ +1 ) , 

where r is the number of pairs eliminated, and ir' G W(k — 2r) is 
maximally crossing in the sense that it cannot be further reduced by 
the procedure of Lemma 1, i.e., for every m we have m </v m + 
Because 7r is crossed, there are at least two pairs left, so k — 2r > 4. 
Consider a pair m ~ 7r / m + t of ir' with minimal £, understood in a 
cyclic sense on the k — 2r positions of n'. Then £ > 2 and all points 
j with m < j < m + £ form pairs with points outside {m, . . . ,m + £}, 
i.e., these pairs cross the pair m ~ 7r / m + £. 

To estimate #S n (ir'), let us first count the choices for indices in 
the interval {m, ...,m + £}, first choosing P rn = (p m ,q m ), P m +e = 
(p m +eiqm+e), and then successively P m+1 = (p m+1 ,q m+1 ), P m+2 , . . .. 
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There are n choices for p m and then, by (CI), o{n 2 ) choices for the three 
indices q m , p m+e , q m+e . Now, p m+1 is fixed by consistency, but q m+ i 
is unconstrained aside from the requirement P m +i T^n P m , which we 
ignore to obtain an upper bound. The same holds for P m +2, ■ ■ ■ , Pm+e-2, 
however, for P m+ £-i there is no freedom due to consistency. For the 
whole interval {m, . . . , m + £} this gives an upper bound n x o{n 2 ) x 
n e ~ 2 = o(n e+1 ) on the number of choices. 

Finally, let us go through the remaining k — 2r — £ + 1 pairs of indices 
in increasing order (cyclically), starting with P m+ e + i. Consistency fixes 
the first index. If m + £+ 1 is paired under tt' with one of the previously 
considered points, then (C2) fixes also the second index within a set 
of size no larger than B. If m + £ + 1 does not form a pair with any 
previous point, then there are at most n choices for the second index. 
As we go through the remaining k — 2r — £ — 1 pairs, the latter case 
occurs exactly — £ times. We deduce that 

#S' n (7r / ) < n r x o(n m ) x n^'* B^' 1 + o(nt +1 ) = o(nt +1 ) . 

□ 

Thus we may restrict the sum in (5) to non-crossing pair partitions of 
{1, . . . , k}, denoted by AfVV(k). For tt G AfVV(k), let PS^n) denote 
the subset of S n (ir) made up of n consistent sequences such that for 
any two sites m,m' paired by tt, m ~„- m', we have P m i = (q m ,Pm), 
where P m = (p m , q m ). Note that for (P 1 , . . . , P k ) e PS n {rr) we have 

(7) E(a n (Pi)-.-a n (P fc )) = l, 

since a n (P m ) = a n {P m i) for equivalent sites m ~ n m' and the a n (P) 
have unit variance. Furthermore, there are not too many sequences 
missed by PS n (Tr): 

Lemma 3. Let tt G MW{k) be non-crossing. Then 

#S n (7T) - #PS n (TT) = Q 

Proof. Let NS n (Tr) = S n (Tr) \ PS n (Tr). Repeating the proof of Lemma 
2, we see that 

#iV,S n (7r) < n x #NS n (Tr') + o(n^ +1 ) , 

where tt' is a partition of k — 2 obtained by removing a nearest neighbor 
pair from tt. Referring to the notation of the proof of Lemma 1, the 
only new feature here is that in case (i), after eliminating the pair 
P m ~n P m +i (and relabeling m + 2 i— > m, etc.), we get a sequence 
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in NS n (n'). This follows because the eliminated pair was of the form 
(P) l) ~« (q,p), so the "defect" is still present in the reduced sequence. 

Iteration of the above gives 

NS n (7r) < J^NS^n") + o(nt +1 ) , 

where n" G VV(2) is a pair partition for which #NS n {ji") = o(n 2 ), by 
condition (C3). □ 

After restricting to pair partitions in (5) and applying Lemmas 1 to 
3 and equation (7), we conclude that 

(8) limEiTr„(l„ fc ) = lim -L £ #P5 n (vr) , 

n^oo fi n—>oo ^-L+o ' ' 

a TreAfVV(k) 

where, by the following lemma, the limit on the r.h.s. is in fact just 
the number of non-crossing pair partitions. 

Lemma 4. If n G AfVV(k), then 

lim #p f" w = 1 . 

712 



n— >0O /n^+l 



Proof. We first note that | — 1 applications of Lemma 1 imply 

#S„(vr) < nf- 1 x #5 n (7r') + o(nf +1 ) , 

where 7r' G W{2) is a pair partition, for which ^S n (n') = n 2 + o(n 2 ), 
by (C3). Thus 

(9) #P5„(tt) < #5 n (7r) < nt +1 + o(nf +1 ). 

In order to prove a lower bound, let us count the number of choices 
for indices, proceeding from left to right. We first have exactly n 2 
choices for This determines p 2 = qi by consistency. If 1 ~„- 2 

then q 2 = pi is fixed since we require a PS n (ir) sequence, otherwise 
there remain exactly n — 1 choices for q 2 (since q 2 = Pi is forbidden 
because 1 </v 2). More generally, once we reach the j th pair of indices, 
Pj = qj-i is fixed by consistency, but qj is fixed if and only if j is paired 
under n with a number i < j. If j is not paired with a prior number, 
there remain at least n — (j — 1)B > n — (k — 1)B choices for qj. 
Indeed this follows from (C2) since qj is free up to the constraint that 
{Pji Qj) l^n (Pi, Qi) for i < j. We obtain in this way the lower bound 

(10) #PS„(7r) > n 2 {n-{k-l)B)"2- 1 > (n - (A; - l)B)^ +l , 
since there are exactly k/2 pairs in tt, and thus k/2 + 1 free indices. 
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Combining the two bounds (9) and (10) completes the proof. □ 

This concludes the proof of Theorem 1, because the number of non- 
crossing partitions, #J\fW(k), is precisely equal to the Catalan num- 
ber Ch [Spe, HP]. ' □ 

2. Examples 

2.1. Flip matrix model. For each n, let $ n : {1, n} — > {1, n} 
be the involution = n+l—p, and take for the equivalence relation 

(P,q) ~n (q,p) ~n ($„(p), ~n ($n(?),$n (?) ) ■ 

Clearly (CI) and (C2) are satisfied. For even n, (C3) holds because 
{Pi q) ~« (q, p') if and only if p = p'. For odd n, this is true except for 
pairs of the form (p, ~ n (^2^> $„(p)), since is a fixed point of 
$ n . As there are only n — 1 such pairs, (C3) holds. Thus, Theorem 1 
applies and the density of states is a semi-circle law. 

Restricting to even n and taking for the matrix elements complex 
Gaussians, we obtain the flip matrix model of [Poi]. Theorem 1 in this 
special case reproduces a result derived in [BMR]. We comment below 
on how this ensemble appears in connection with the Anderson model. 

2.2. Condition (C3). The importance of condition (C3) is illustrated 
by the following example. Let the equivalence relation be induced by 

(P,q) ~n (q,p) ~n ($n(p),9) ~n ($ n {q),$ n (p) ) , 

where <3> n is as above. Unless p or q are fixed points of $ n , the class 
of (p, q) contains 8 elements. Conditions (CI) and (C2) are clearly 
satisfied. However, (C3) is violated since (p,q) ~ n (q,®n(p)) f° r every 
p,q, giving at least n x (n — 1) pairs (p, q) ~ n (q,p') with p ^ p'. 
Consider, associated to this ~ n , an ensemble of real symmetric matrices 
with a n (P) = a n (P') when P ~ n P' . The matrix X n has a [n/2] 
dimensional null space (here [a] denotes the integer part of a e M), 
consisting of vectors (t>i, . . . , i> n ) with Vj = —v n+ i-j. As a result, the 
density of states for this ensemble has an atom at 0. In fact, the density 
of states can be computed by noting that 

where = denotes unitary equivalence and F n is an ensemble of matrices 
satisfying the conditions of Theorem I. Thus the density of states is a 
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superposition of an atom at and a scaled semi-circle law: 

- S(x) dx + - — V4 - 2a; 2 X {x 2 < 2) V2 dx . 
2 2 2ir 



2.3. Higher dimensional flip models. Let us first describe how 
these models appear as effective Hamiltonians for a finite size approx- 
imation of the Anderson model in the weak coupling limit, and then 
apply Theorem 1 to them. 

The Anderson Hamiltonian H = H + V on the £ 2 space over a square 
lattice is the sum of a translation invariant part H = Hq and a (small) 
random potential V. Consider such an operator in a finite volume 

f L L 

= {-2 +1 "-"2 

together with periodic boundary conditions. Then V is a multiplica- 
tion operator on £ 2 (Al) given by a random real- valued function v on 
Al- Furthermore the orthonormal basis <j> p (x) = e ip ' x / L d / 2 , with quasi- 
momentum p in the (discrete) Brillouin zone = ^-A L , consists of 
eigenfunctions of H with eigenvalues £{p), where S is some smooth 
periodic function on the torus [— tt, n) d — the so-called symbol of H . 

For small V, the spectral analysis of H near a fixed energy E involves 
primarily states P with quasimomenta in the discrete Fermi surface 

(11) A2(£) = {peA* L \\£(p)-E\<^}, 

with f3 > a parameter giving its width. The restriction of the random 
potential to the corresponding subspace will be the higher dimensional 
flip model. It is hence the lowest order in a nearly degenerate pertur- 
bation theory for the Anderson Hamiltonian at scale L. 

More precisely, the multiplication operator ip(x) h- > v(x)ip(x) acts on 
the basis P as follows: 

<j)p ^ - Q)<i>q » 

where v is the (finite sum) Fourier transform of v, i.e., 

x&A L 
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and p — q is evaluated modulo 2tt. Restriction to the Fermi surface 
gives thus the following convolution operator A w on £ 2 (A* L (E)) 

A u (f>(p) = Yl v(p - q)(j>{q) ■ 
q eA* L (E) 

Hence we are led to introduce the random matrix 

(12) Y n (p,q) = v(p-q) , p,q e A* L (E) . 

Here n = #A* L (E) = 0(L d - r ). 

If one takes v(x) to be independent and identically distributed cen- 
tered Gaussians, then the variables v{p) are also centered Gaussians, 
which are mutually independent apart from the constraint v(p) = 
v(—p) which results because the v(x) are real. In terms of the variance 
A 2 = E(v(x) 2 ), the variance of v(p) is X 2 /L d . The higher dimensional 
flip matrix model is the random matrix Y n defined in (12), with the 
index set given by (11), where random variables v{p) have variance 
\ 2 /L d and satisfy the constraint v(p) = v(—p), as well as the bounds 

max E(|%))| fe ) < M k \ k L~^ . 

P eA* E 

To simplify the discussion, we suppose that £(p) = p 2 /2 in a neighbor- 
hood of E > 0, so that A* L (E) is the intersection of a spherical shell 
with the lattice A* L . This eases some arguments, but does not affect 
the end result provided the level sets {p E A* L | £(p) = E±5}, for small 
S, are sufficiently regular and have curvature bounded away from zero. 

Remark: The flip matrices of Section 2.1 are a simplification of Y n in 
d = 2 obtained via a course graining procedure. The shell A* L (E) is 
divided into m blocks, each subtending an angle of size 0{L-2). Each 
block contains 0(L~2 x L~ L x L 2 ) = 0(L2) points. The flip matrix 
F m (b,b f ) is indexed by pairs of blocks, with Gaussian E m (b,b') which 
are pairwise independent except that 

F m (b,b') = F m {b',b) = F m (-b',-b) = F m (-b, -b>) . 

The above coarse graining was introduced in [Poi, BMR] as a simplifi- 
cation to avoid the combinatorics required to deal with matrix elements 
Y n (p, q) with p — q small. We do not need coarse graining here, because 
Theorem 1 applies directly to an appropriate rescaling of Y n . 

We now turn to the analysis of the higher dimensional flip model. 
The matrix elements Y n (p, q) and Y n (p', q') are statistically independent 
unless p — q = ±(p' — q'), and thus we define ~ n accordingly: 

(p, q) ~n (p 1 , q') p-q = ± (p' - q') 
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Condition (C2) is clearly satisfied, but conditions (CI) and (C3) require 
proof. 

Condition (C3) is a bit easier to verify. The relevant pairs (p, q) ~ n 
(q,p') with p 7^ p' have to satisfy p' = 2q — p. Fixing q G A* L (E), this 
implies that the conditions p G A* L (E) and 2q—p G A* L (E) have to hold 
simultaneously. The relevant set of p's is contained in the intersection 
of two spherical shells: 

\p 2 /2-E\ < f \ n [\{p-2q) 2 /2-E\ < | 



Lj {'^ ^' ' L 

As \q\' 2 /2 = E, these shells are barely touching and there are 

0(L-^L d ) = O(L^) = o(L d - 1 ) = o(n) 

points in the intersection of this set with A* L (E), essentially because the 
(i-dimensional volume of the set is O^L^^) and the lattice cell volume 
of is (2ir) d / L d . As this holds for all n values of q, we conclude that 
(C3) holds. 

For (CI), we need to estimate, given a fixed value of p, the number 
of triples q,p',q' G A* L (E) with (p, q) ~ n {p',q'). First note that 

q' = p' ± (p-q) , 

so that q' takes 2 values once q and p' are chosen. Furthermore 

\qf-\pf = ±2p'-(p-q) + b-g| 2 . 

Because p',q' G A* L (E), this implies 

2(3 



\p-q\ , / (p-?) 
— r — ± p 



\p 



< 



L\p-q\ 



In order to count the number of possible values for q and p', we 
consider separately the cases of large and small p — q. If \p — q\ > L~2 , 
then 



\p-q\ 



2|p'| 



± cos^ 



< 2 " 



Ip'IvT 



where 8 is the angle between p' and {p — q). Because for all t 
\{9 G [-7r,7r] I | cos 6* - t\ < 5}\ < V25 , 

we conclude that the allowed set of p' is contained in a spherical 
cone which subtends an angle no larger than 0(L~*). The intersec- 
tion of this spherical cone with A* L (E) contains O^- 1 ^- 1 ) x L d ) = 
0(L-^ d ~ 1] ) points. As there are 0(L d - v ) points q G A* L (E) with 
\p — q\ > L~2 , we end up with no more than 0{Li^ d ^) triples in this 
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case. However, in the other case there are only OiL- 1 —^- 1 ^ x L d ) = 
0(L^ d ~^) points q G A* L (E) with \p — q\ < L~~z, and thus no more 
than 0(L^ d -^ x L d ~ r ) = 0(L^ d -^) triples. Altogether we conclude 
that 

and thus that (CI) holds. 

Therefore, Theorem 1 applies to X n = c n Y n with a scalar c n chosen 
so that the variance of the matrix elements of X n is 1/n, namely c 2 = 
(L d /\ 2 ) x n~ x = 0(L)/\ 2 . In the context of the Anderson model, it is 
natural [Poi, Bel] to take L A~ 2 , giving n(A) = #A* L (E) = 0{\ 2 ~ 2d ). 
Then Y n ^ (p, q) A 2 x X n ^\) ip, q) where the X n are matrices which 
satisfy the hypotheses of Theorem 1. It follows that the density of 
states of Y n (x) is asymptotic as A — > to a scaled semi-circle law of 
width of order A 2 . 



3. Asymptotic Freeness 

The semicircle law is the only stable law of Voiculescu's free proba- 
bility theory for which all moments exist, playing a role there analogous 
to that of the Gaussian in classical probability. The above techniques 
also allow to extend known results on asymptotic freeness of random 
matrices [VDN, Voi, Spe, HP]. Suppose given a family (X$) ieI of ran- 
dom n x n matrices with some index set I. For different indices i, the 
entries of the matrices are supposed to be independent, and for each 
i, the family (Xn^)neN is specified by a possibly i-dependent equiva- 
lence relation satisfying conditions (CI), (C2) and (C3). Furthermore 
suppose given a second family (Dn ^)j e j of diagonal real random ma- 
trices, the entries of which are also independent among themselves and 
with respect to the family (X$)iei- Each (_D^) ngN is supposed to be 
such that a well-defined density of states exists, namely a probability 
measure A/} such that 

lim E — Tr n (f(D^)) = [ Afj(dx) f(x) , / G C (K) . 

n— >oo n J 
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The asymptotic behavior (in n) of the families (Xn^)iei and (Dn^)jeJ 
now defines a state <p on the algebra A = C[{x®} ieI , {d {j) } j&J \ of non- 
commutative polynomials by 

= lim E — Tr n (Q({X^} ieI , {D®} jeJ )) , QeA. 

n^oo n 

Theorem 2. The commutative subalgebras C[x^] anc? C[d^] are free 
in (A, (fi) for all i G I and j G J. Moreover, the marginals of <p on 
C[x^] and C[d^] are respectively the semicircle law and the measure 
A/}. 

The first statement means, in particular, that all mixed moments 
of tp can be calculated from the moments of the semicircle law and 
the densities A/} by use of free combinatorics discussed in detail in 
[Spe, HP]. In view of the above arguments, the proof of this theorem 
is a straightforward generalization of the one given in [Spe]. 

Acknowledgment: We would like to thank J. Bellissard, for bringing 
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